Search results for "Adaptive mesh refinement"
showing 10 items of 23 documents
A study on time discretization and adaptive mesh refinement methods for the simulation of cancer invasion: The urokinase model
2016
In the present work we investigate a model that describes the chemotactically and proteolytically driven tissue invasion by cancer cells. The model is a system of advection-reaction-diffusion equations that takes into account the role of the serine protease urokinase-type plasminogen activator. The analytical and numerical study of such a system constitutes a challenge due to the merging, emerging, and traveling concentrations that the solutions exhibit. Classical numerical methods applied to this system necessitate very fine discretization grids to resolve these dynamics in an accurate way. To reduce the computational cost without sacrificing the accuracy of the solution, we apply adaptive…
Simulations of stellar/pulsar wind interaction along one full orbit
2012
The winds from a non-accreting pulsar and a massive star in a binary system collide forming a bow-shaped shock structure. The Coriolis force induced by orbital motion deflects the shocked flows, strongly affecting their dynamics. We study the evolution of the shocked stellar and pulsar winds on scales in which the orbital motion is important. Potential sites of non-thermal activity are investigated. Relativistic hydrodynamical simulations in two dimensions, performed with the code PLUTO and using the adaptive mesh refinement technique, are used to model interacting stellar and pulsar winds on scales ~80 times the distance between the stars. The hydrodynamical results suggest the suitable lo…
Adaptive discontinuous evolution Galerkin method for dry atmospheric flow
2014
We present a new adaptive genuinely multidimensional method within the framework of the discontinuous Galerkin method. The discontinuous evolution Galerkin (DEG) method couples a discontinuous Galerkin formulation with approximate evolution operators. The latter are constructed using the bicharacteristics of multidimensional hyperbolic systems, such that all of the infinitely many directions of wave propagation are considered explicitly. In order to take into account multiscale phenomena that typically appear in atmospheric flows nonlinear fluxes are split into a linear part governing the acoustic and gravitational waves and a nonlinear part that models advection. Time integration is realiz…
Spectral WENO schemes with Adaptive Mesh Refinement for models of polydisperse sedimentation
2012
The sedimentation of a polydisperse suspension with particles belonging to N size classes (species) can be described by a system of N nonlinear, strongly coupled scalar first-order conservation laws. Its solutions usually exhibit kinematic shocks separating areas of different composition. Based on the so-called secular equation [J. Anderson, Lin. Alg. Appl. 246, 49–70 (1996)], which provides access to the spectral decomposition of the Jacobian of the flux vector for this class of models, Burger et al. [J. Comput. Phys. 230, 2322–2344 (2011)] proposed a spectral weighted essentially non-oscillatory (WENO) scheme for the numerical solution of the model. It is demonstrated that the efficiency …
Fine-Mesh Numerical Simulations for 2D Riemann Problems with a Multilevel Scheme
2001
The numerical simulation of physical problems modeled by systems of conservation laws can be difficult due to the occurrence of discontinuities and other non-smooth features in the solution.
Adaptive mesh refinement techniques for high-order shock capturing schemes for multi-dimensional hydrodynamic simulations
2006
The numerical simulation of physical phenomena represented by non-linear hyperbolic systems of conservation laws presents specific difficulties mainly due to the presence of discontinuities in the solution. State of the art methods for the solution of such equations involve high resolution shock capturing schemes, which are able to produce sharp profiles at the discontinuities and high accuracy in smooth regions, together with some kind of grid adaption, which reduces the computational cost by using finer grids near the discontinuities and coarser grids in smooth regions. The combination of both techniques presents intrinsic numerical and programming difficulties. In this work we present a …
An adaptive rectangular mesh administration and refinement technique with application in cancer invasion models
2022
We present an administration technique for the bookkeeping of adaptive mesh refinement on (hyper-)rectangular meshes. Our technique is a unified approach for h-refinement on 1-, 2- and 3D domains, which is easy to use and avoids traversing the connectivity graph of the ancestry of mesh cells. Due to the employed rectangular mesh structure, the identification of the siblings and the neighbouring cells is greatly simplified. The administration technique is particularly designed for smooth meshes, where the smoothness is dynamically used in the matrix operations. It has a small memory footprint that makes it affordable for a wide range of mesh resolutions over a large class of problems. We pre…
Error Estimates and Automatic Adaptive Mesh Refinement for the Metal Forming FEM Analysis
1988
The Authors propose a new technique which enables a estimation of the error inherent with the FEM analysis of metal forming processes. The aim is to evaluate the zones where the error is higher in order to proceed to a refinement of the mesh in such zones, and to obtain a smaller value of the global error. Moreover, to simplify the analyst work in the progressive refinement of the mesh, it has been prepared a software able to read the drawing created by a CAD program and to generate, automatically, all the geometrical and topological data necessary to perform the analysis on Personal Computer. The automatic renumbering of the elements in the refined mesh has been performed with the aim to r…
Comparison between adaptive and uniform discontinuous Galerkin simulations in dry 2D bubble experiments
2013
Accepted by the Journal of Computational Physics Adaptive mesh refinement generally aims to increase computational efficiency without compromising the accuracy of the numerical solution. However it is an open question in which regions the spatial resolution can actually be coarsened without affecting the accuracy of the result. This question is investigated for a specific example of dry atmospheric convection, namely the simulation of warm air bubbles. For this purpose a novel numerical model is developed that is tailored towards this specific application. The compressible Euler equations are solved with a Discontinuous Galerkin method. Time integration is done with an IMEXmethod and the dy…
A semi-Lagrangian AMR scheme for 2D transport problems in conservation form
2013
In this paper, we construct a semi-Lagrangian (SL) Adaptive-Mesh-Refinement (AMR) solver for 1D and 2D transport problems in conservation form. First, we describe the a-la-Harten AMR framework: the adaptation process selects a hierarchical set of grids with different resolutions depending on the features of the integrand function, using as criteria the point value prediction via interpolation from coarser meshes, and the appearance of large gradients. We integrate in time by reconstructing at the feet of the characteristics through the Point-Value Weighted Essentially Non-Oscillatory (PV-WENO) interpolator. We propose, then, an extension to the 2D setting by making the time integration dime…